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Introduction 


The  goal  of  this  research  project  is  to  investigate  the  relationship  between  primary  breast  tumor 
vascular  development  and  metastatic  potential.  To  facilitate  this  research  in  an  in  vivo  mouse 
model,  we  have  developed  a  novel  dual  modality  system  for  small  animal  imaging.  A  high 
sensitivity  optical  imaging  component  is  used  to  detect  bioluminescent  light  emitted  from  implanted 
tumor  cells  that  have  previous  been  transfected  with  the  luciferase  gene.  In  the  same  instrument, 
these  high-sensitivity  cameras  can  be  coupled  with  scintillators  and  collimators  to  allow  Single 
Photon  Emission  Computed  Tomography  imaging.  While  the  bioluminescent  component  will  be 
sensitive  to  the  detection  of  distant  metastases,  the  nuclear  medicine  aspect  of  the  system  will  be 
utilized  in  subsequent  years  for  studying  tumor  angiogenesis  and  vasculature. 

Body 

Our  statement  of  work  specified  that  the  first  year  would  be  dedicated  to  the  development  of  a 
fully  integrated  multi-head  optical/BLI  and  micro-SPECT  system.  Task  1  was  subdivided  into  8 
steps,  1-h.  In  summary: 

a-  Gantry  and  camera  housings  have  been  designed,  built  and  tested  . 

b-  Camera  and  support  motion  assemblies  were  built. 

c-  The  CCDs  and  associated  first  generation  optical/scintillator  components  were  built  and 
tested. 

d-  The  data  acquisition  system  has  been  completed  and  assembled  with  the  above  components. 
These  have  been  used  with  a  test  optical  imaging  system  based  on  SITe  CCD  chips,  now 
completed  and  in  use. 

e-  The  5-head  system  was  constructed,  but  we  experienced  CCD  chip  failure.  As  explained 
below,  this  caused  us  to  move  to  a  new  CCD  assembly,  with  the  same  price  but  different 
characteristics.  This  required  a  redesign  of  the  electronics,  still  in  progress. 

f-  Single  camera  tests  were  successful,  showing  high  sensitivity. 

g-  Development  of  optical  MLEM  reconstruction  software  has  progressed,  and  we  have 
reported  our  first  results  in  tomographic  imaging  of  light  emission.  Our  images  show  that  it 
is  possible  to  detect  small  lesions  not  only  subcutaneously  but  also  at  depth  in  vivo,  and  we 
characterizef  the  advantages  of  this  method  for  quantitative  imaging.  3D  images  were 
obtained  and  reported  at  the  ERA  of  HOPE  symposium. 

h-  Phantom  testing  has  been  completed  including  both  homogeneous,  cylindrical,  and 
heterogeneous,  mouse-shaped  phantoms.  We  have  also  demonstrated  micro-SPECT  (Single 
Photon  Emission  Computed  Tomography)  imaging  at  high  resolution. 

We  report  below  on  our  progress  along  these  points.  While  the  completion  of  an  integrated 
optical/microSPECT  imaging  device  was  temporarily  hampered  by  the  commercial  failure  of  the 
SITe  manufacturer,  we  selected  new  CCDs  with  equivalent  or  better  performance.  The  completion 
of  the  new  electronics  is  now  advanced.  Because  of  the  new  features,  we  expect  that  the  new  device 
will  result  in  portable  devices  suitable  for  human  applications  at  a  later  time. 
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High  Sensitivity  CCD  camera  (a-f) 


The  dual  optical/SPECT  imaging  system  uses  CCD  (charge-coupled  device)  cameras 
selected  for  their  high  and  nearly  constant  sensitivity  over  the  full  range  of  wavelengths  commonly 
used  in  optical  imaging,  from  blue  to  near  infrared.  Based  on  our  design  the  required  parameters  of 
the  CCD  are:  non-color,  back-illuminated,  full  frame  image  sensor,  512x512  pixels,  with  large  pixel 
size  to  provide  good  dynamic  range.  Preferably  the  quantum  efficiency  of  the  CCD  should  be 
greater  than  85%  from  400  nm  to  750  nm,  and  should  remain  above  50%  up  to  900  nm.  The 
originally  selected  CCD  SI-032AB  produced  by  Scientific  Imaging  Technologies  Inc.,  Tigard,  OR 
fulfilled  all  this  requirements  but  it  is  no  longer  available  due  to  the  company’s  bankruptcy. 

Thus,  a  new  CCD  was  selected,  CCD77-00  (e2v  technologies,  inc.  4  Westchester  Plaza,  PO 
Box  1482,  Elmsford  NY  10523- 1482  USA).  This  CCD  has  a  pixel  size  of  24x24  pm  providing  a 
large  well  capacity  of  300,000  e",  with  a  sensitivity  of  2.5  pV/e~,  low  dark  current  (20  pA/cm2  at  20 
°C),  and  low  readout  noise  (3  e  RMS)  providing  a  dynamic  range  of  ~1 00,000.  The  CCD  is  cooled 
to  -40°  C  reducing  the  dark  current  signal  to  less  than  0.1  eVpixel/sec.  All  relevant  parameters  are 
similar  or  better  than  those  of  SI-032AB,  but  the  support  electronics  requires  significant  changes. 
We  are  currently  implementing  the  new  driving  circuits  (Fig.l)  which  can  provide  the  new  control 
sequence  and  voltages. 


Figure  1.  The  implemented  electronic  circuitry  for  driving  the  CCD  signals. 

The  CCD  is  incorporated  in  a  self-contained,  cooled  camera,  equipped  with  electronic  circuitry  and 
fast  optics  (25  mm  focal  length,  170.95).  Each  camera  is  calibrated  using  a  low-intensity,  diffuse,  flat 
field  source  that  gives  a  known  radiance  (typically  3.0  x  10  7  W/cm2/sr).  The  light  source  is 
periodically  checked  for  uniformity  using  a  NIST-traceable  research  radiometer  (model  IL  1700, 
International  Light,  Inc.  Newburyport,  MA).  By  imaging  this  source  the  digital  units  provided  by 
the  camera  digitizer  can  be  converted  directly  into  absolute  physical  units  (W/cm2/sr  or 
photons/sec/cm2/sr).  This  method  accounts  for  the  transmission  efficiency  of  the  entire  optical 
system,  and  corrects  the  non-uniformity  in  the  field  of  view  due  to  lens  vignetting  and  variations  in 
pixel  sensitivity. 

Multiple  Head  Optical  Imaging  System 
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Multiple  high  sensitivity  CCD  cameras  simultaneously  record  views,  and  a  computer  controlled 
rotation  mechanism,  allows  imaging  at  multiple  angular  positions,  as  required  for  three  dimensional 
reconstructions.  The  support  electronics  allow  simultaneous  control  of  the  cameras  for  light 
exposure,  image  readout  and  preprocessing,  and  temperature  and  vacuum  control.  A  horizontal  bed, 
made  of  sparse  mesh  material  to  reduce  interference  with  light,  is  used  to  support  and  immobilize 
the  animal  during  imaging.  A  gas  anesthesia  unit  is  connected  to  the  bed  and  to  the  animal.  To 
exclude  ambient  light,  the  system  is  encapsulated  in  a  light-tight  enclosure.  Light  images  at  each 
position  of  the  gantry  are  acquired  for  co-registration  with  the  bioluminescence  image  using  a  set  of 
6  diffuse  light  sources. 


Figure  2.  The  design  of  the  imaging  system  includes  5  high  sensitivity  cameras  mounted  on  a 
rotating  gantry,  and  a  transparent  animal  positioning  system. 


Validation  studies  using  phantoms  verified  the  utility  of  the  approach  and  accuracy  of  the  3D 
reconstruction  algorithms.  Here  we  describe  only  two  phantoms  of  increasing  complexity,  a 
homogeneous  diffusing  cylinder  and  a  homogeneous  mouse-shaped  phantom. 

Homogeneous  cylindrical  phantom  (h) 

Figure  3a  shows  the  experimental  setup  with  a  30  mm  diameter  cylindrical  phantom  filled  with  1  % 
Intralipid  /  1%  Agarose  gel  (Intralipid  gel  for  brevity),  which  approximates  light  scattering  in 
tissues.  An  optical  fiber  (1  mm  diameter)  was  placed  in  the  phantom  before  congealing  the  filling 
mixture  and  was  optically  coupled  to  a  560nm  LED.  A  set  of  20  images  was  obtained  for  each 
experiment  by  rotating  the  phantom  18°  at  a  time. 

The  total  amount  of  light  captured  at  each  imaging  angle  both  in  air  and  Intralipid  phantom  is 
presented  in  Figure  3b,  obtained  while  keeping  the  light  source  intensity  constant.  The  fiber  optic 
source  in  air  shows  a  small  variation  due  to  changes  in  distance  from  the  lens.  When  the  source  was 
immersed  in  Intralipid  it  was  noted  that  for  angles  between  0  and  180°,  when  the  light  source  is 
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closest  to  the  camera,  up  to  4  times  more  light  was  captured  than  for  the  source  in  air,  as  the 
scattering  of  light  in  the  phantom  caused  a  significant  increase  in  the  photon  flux  toward  the 
objective.  When  the  source  was  in  the  position  furthest  from  the  camera  (270°)  only  half  as  much 
light  was  captured  because  the  light  source  was  shadowed  by  the  body  of  the  phantom.  Overall,  a 
variability  of  almost  an  order  of  magnitude  was  observed  when  imaging  this  phantom  at  different 
angles,  reducing  the  precision  and  reliability  of  a  quantitative  assessment  from  single  planar  images. 
This  effect  would  be  expected  to  affect  animal  studies,  when  light  is  emitted  at  depth  within  the 
body  from  multiple  sources,  such  as  small  metastatic  tumors  in  the  lungs  or  other  internal  organs. 

Figures  3c  and  3d  show  images  obtained  first  in  air  and  then  in  gel,  with  the  fiber  displaced 
10  mm  from  the  center  at  a  viewing  angle  of  180°.  A  cross-section  through  the  reconstructed  3D 
image  is  also  shown  for  each  case.  The  position  of  the  source  was  correctly  identified  by  the 
reconstruction  algorithm.  The  total  light  signal  in  air  and  Intralipid  differed  by  only  3%  (1,435  x 
107  photons/sec  for  air  vs.  1,390  x  107  photons/sec  for  Intralipid),  further  showing  the  quantitative 
nature  of  the  3D  reconstruction  algorithm  in  the  presence  of  diffusion.  Confirmatory  results  were 
also  obtained  for  multiple  light  sources,  both  in  air  and  Intralipid,  e.g.,  for  two  sources  separated  in 
depth  by  10  mm  (Fig.  3e). 


Figure  3.  Validation  of  LET  in  phantoms. 

a.  Experimental  setup:  30  mm  cylindrical  phantom  filled  with  1%  Intralipid  gel. 

b.  Total  amount  of  light  captured  from  a  fiber  optic  source  in  air  and  in  the  phantom 
at  various  angles.  The  light  source  is  displaced  10  mm  from  the  axis  of  rotation. 


7 


c,  d.  Image  of  light  emitting  fiber  in  air  (c)  and  in  intralipid-filled  medium  (d)  displaced  10  mm  from  the 
center.  The  intensity  profiles  through  the  reconstructed  3D  image  are  also  shown. 

e.  3D  intensity  profiles  for  1%  Intralipid  phantom  with  two  sources  separated  by  10  mm. 

Homogeneous  mouse  shaped  phantom 

In  order  to  investigate  the  LET  reconstruction  in  the  context  of  a  complex  body  surface,  a  mouse 
shaped  1%  Intralipid  gel  phantom  with  a  9.5  mm  diffuse  spherical  light  source  embedded  in  the 
thoracic  region  was  completed  and  imaged  on  the  LET  system.  Figure  4  shows  a  longitudinal 
section  of  the  phantom  with  the  3D  reconstructed  source  overlay.  The  reconstructed  source  is  10 
mm  in  size,  in  agreement  with  its  physical  dimensions,  and  at  the  correct  location  (within  1mm). 


Figure  4.  Longitudinal  section  of  the  mouse  shaped  phantom  with  the  3D  reconstructed  source 

overlay. 


Optical  MLEM  Software  (f) 

Three-Dimensional  Bioluminescent  Imaging  (3D-BLI)  of  luciferase-expressing  cells  in  live 
small  animals  would  greatly  enhance  applications  in  biomedicine  since  light  emitting  cell 
populations  could  be  unambiguously  associated  with  specific  organs  or  tissues.  We  consider  this 
pivotal  for  investigating  tumor  growth,  metastasis,  and  specific  biological  molecular  events.  The 
imaging  approach  must  account  for  the  main  optical  properties  of  biological  tissue  because  light 
emission  from  a  distribution  of  sources  at  depth  is  strongly  attenuated  due  to  optical  absorption  and 
scattering  in  tissue.  Our  image  reconstruction  method  is  based  on  the  deblurring  Expectation 
Maximization  (EM)  method  and  takes  into  account  both  of  these  effects.  To  determinate  the 
boundary  of  the  object  (skin  surface),  we  use  the  standard  iterative  algorithm  -  Maximum 
Likelihood  (ML)  reconstruction  method.  Depth-dependent  corrections  for  quantitative  light 
intensity  were  included  in  the  reconstruction  procedure  by  using  the  diffusion  equation  for  light 
transport  in  semi-infinite  turbid  media  with  extrapolated  boundary  conditions. 

For  detection  of  bioluminescence  emitting  probes  in  living  animals,  it  is  necessary  to  use 
high  sensitivity  low-noise  detectors,  such  as  cooled  Charged-Coupled  Devices  (CCDs),  and 
advanced  software  codes  for  image  reconstruction.  Reconstruction  algorithms  must  take  into 
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account  the  main  optical  properties  of  biological  tissue:  light  absorption,  scattering,  and  possibly 
reflection  at  major  interfaces.  In  this  case,  a  model  describing  light  transport  in  turbid  media  is 
critical  for  any  reconstruction  technique15 .  Widely  used  methods  include  the  Monte  Carlo 
approach10-19  and  the  Radiative  Transfer  Equation  (RTE)  20-24 . 

The  Monte  Carlo  method  is  a  discrete  model  of  individual  photon  interactions  that  can  be 
described  by  computer  simulation  of  appropriately  weighted  random  absorption  and  scattering 
events.  This  method  is  an  accurate  and  powerful  way  to  study  photon  transport  in  tissue  requiring 
no  simplifications,  but  expensive  computational  resources  are  often  needed  due  to  the  slow 
statistical  nature  of  the  method. 

The  RTE  method  is  a  continuous  model  based  on  an  integro-differential  equation, 
specifically  a  Boltzman-like  transport  equation  for  the  photon  wave  intensity,  which  includes 
emission,  absorption,  and  scattering  in  the  media.  This  equation  most  often  can  not  be  solved 
analytically,  and  is  often  simplified  by  a  diffusion  approximation 20  for  numerical  solution.  The 
diffusion  equation  gives  us  not  only  a  great  simplification  of  the  problem  but  a  practical  tool  to 
describe  the  diffuse  part  of  the  radiative  intensity.  A  disadvantage  of  this  method  is  the  inaccuracy 
in  predicting  the  light  distribution  near  light  sources  and  boundaries.  The  diffusion  approximation  is 
only  valid  for  materials  where  scattering  dominates  over  absorption. 

FORWARD  MODEL:  LIGHT  DIFFUSION  IN  TURBID  MEDIA 

To  determine  the  position  of  bioluminescence  sources  we  must  describe  photon  propagation 
through  tissue  (light  can  be  considered  as  a  stream  of  photons)  and  simplify  the  RTE  to  the 
diffusion  approximation.  The  RTE  is  an  integro-differential  equation  derived  by  considering  energy 
conservation  in  a  small  volume  dV.  The  photon  distribution  function  N(r,x  ,  t)  depends  on  time  t 
and  can  be  spatially-varying  with  any  position  r  (x,  y,  z)  in  direction  s(6,<p).  L(r, s ,  t)  =  c 
h  v  N(r,  s ,  t)  is  the  radiance  -  power  per  unit  area  per  unit  solid  angle  Q ,  h  is  the  Planck  constant,  v 
is  photon  frequency,  and  c  is  photon  speed  in  vacuum.  Photons  can  only  be  added  or  subtracted 
from  the  photon  distribution  function  if  they  interact  in  this  volume  dV  through  the  following 
processes20 :  i)  photons  are  lost  through  the  boundary  of  media  with  different  refraction  index  n ;  ii) 
photons  incoming  in  direction  s  can  be  absorbed  with  absorption  coefficient  jua  (r)  or  scattered  to 

any  another  direction  s'  with  scattering  coefficient  pt  (r)  and  a  probability  (phase  function) 
p(s,s') ;  iii)  photons  can  be  gained  through  scattering  from  any  direction  s'  into  the  direction  x; 
iv)  photons  are  also  gained  through  an  internal  light  source  S(r ,s,  t). 

If  we  neglect  interference  effects  of  photons  and  drop  the  integrals  over  volume,  we  can  write 
the  time-dependent  RTE  for  the  radiance: 


(1) 


where  cM  =c/n  is  the  speed  of  light  in  the  tissue. 

By  applying  the  PI  expansion  approximation  and  assuming  that  the  scattering  coefficient 
dominates  the  absorption  coefficient,  tire  RTE  can  be  further  simplified  to  a 
time-dependent  diffusion  equation: 
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(4) 


t 


—  =  -OV20(r,  t)  -  Ma<p{r,  t)  +  S(r,  t) 

cn  dt 

The  analytic  solution  for  an  isotropic  point-like  source  S(r,t )  =  5{r,t )  is  the  time-domain 
Green's  function: 

G(r,  t )  =  cn  (4  nDcnyiJ2  (5) 

and  the  solution  to  the  diffusion  equation  (4)  for  radiance  is  given  by  the  convolution: 

(p{r,t)  =  G(r,t)®8(r,t)  (6) 

For  steady-state  diffusion,  as  is  the  case  with  BLI,  we  can  neglect  the  time-dependence  in 
equation  (4)  and  obtain: 

V(r)  =  V'"'  (7) 

r 

where  the  source  power  S0  is  constant  in  amplitude,  and  the  radiance  decreases  exponentially 
with  distance  r,  according  to  an  effective  attenuation  coefficient//^  =  (//a  ID)'12 . 

The  next  problem  is  the  refraction  of  photons  at  the  tissue-air  interface  that  will  reduce  the 
fluence  rate  out  of  the  animal.  We  have  adopted  the  solution  for  semi-infinite  medium  with 
extrapolated  boundary  and  photon  dipole  source  model 2  3 1-34  and  used  it  to  quantify  the  light 
intensity  from  experimental  data  for  phantom  filled  with  Intralipid j5 .  This  diffusion  solution  has 
been  shown  to  be  in  good  agreement  with  Monte  Carlo  simulations  and  experiments2  31,34 . 

The  diffusion  approximation  is  very  useful  for  estimating  the  surface  photon  intensity  and 
spot  size  depth-dependence  for  BLI  experiments,  especially  for  in  vivo  small  animal  studies.  Figure 
5  shows  surface  radiance  and  its  spot  size  as  a  function  of  depth  for  different  tissue  optical 
parameters.  The  distributions  for  the  absorption  surface  intensity  for  several  sets  of  absorption 
coefficients  are  plotted  in  Fig. la,  while  the  surface  intensity  for  variable  scattering  coefficients  is 
plotted  in  Fig. lb.  Figurelc  shows  the  parameter  Full  Width  Half  Maximum  (spot  size)  for  different 
concentration  of  Intralipid:  higher  concentrations  of  Intralipid  give  lower  photon  intensity  level,  but 
spatial  resolution  is  better  for  depths  -1-1.5  cm.  This  depth  is  consistent  for  good  results  for  small 
animal  BLI  reconstruction. 


depth.  cm 


depth,  an 


depth,  or 


a  be 

F'g-5  Surface  Intensity  Peak  (a)-(b)  and  FWHM  (c)  as  a  function  of  source  depth 
calculated  for  different  concentration  of  Intralipid  by  diffusion  approximation  with 
extrapolated  boundary  condition. 
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Now  we  can  calculate  distributions  of  surface  intensity  for  two  sources  with  some  separation 
as  a  function  of  source  depth.  As  shown  in  Fig.6  intensity  decreases  with  depth,  and  at  more  than  12 
mm  depth  the  two  peaks  are  practically  inseparable.  In  this  case,  reconstruction  of  two  sources 
becomes  problematic. 


X.  nvn  X,  wn  X,  im 


a  b  c 

Fig.6  Surface  Intensity  for  two  sources  with  separation  10  nun  as  a  function  of  depth:  a)  4  mm, 
b)  8  mm  and  c)  12  mm  for  1%  Intralipid  concentration. 

IMAGE  RECONSTRUCTION  METHODS 


We  have  developed  a  3-D  reconstruction  algorithm  for  small  animal  imaging,  and  this  paper 
presents  the  results  for  reconstruction  in  the  multi-scattering,  continuous  wave  regime.  By  using 
detectors  with  constant  response  over  the  entire  range  of  wavelengths  expected  for  the  light- 
emission  problem  at  hand,  our  images  are  proportional  to  the  number  of  quanta  captured  in  each 
pixel.  Our  reconstruction  approach  consists  of  two  steps:  light  surface  reconstruction  to  determine 
the  geometry  of  turbid  media  boundary  followed  by  source  reconstruction  inside  the  tissue. 

The  surface  reconstruction  is  based  on  the  ML  method 40-4 ' .  Using  the  ML  algorithm,  a 

til 

sequence  of  image  estimates  is  generated  for  the  k  iteration: 
nT  =  where  =Tuavnkj 


(12) 


£+1  k 

where  rij  and  rij  are  image  intensity  values  in  the  new  and  current  estimate  for  the  voxel  j  . 

The  measured  data  mj  collected  along  lines  i  which  connect  a  voxel  on  the  object’s  surface  to  the 

k 

CCD  camera,  gave  values  proportional  to  the  light  intensity  detected  by  the  pixel;  q .  is  the  expected 

k 

count  in  line  /  for  the  current  estimate  n  ,  and  a.,  is  the  probability  that  a  light  quantum  emitted  at 

V 

t 

voxel  j  will  be  detected  in  the  line  i.  The  normalization  factor  ^joij  includes  all  possible  measured 

/=i 

lines.  Figure  7  shows  examples  of  ML  surface  reconstruction  for  the  phantom  with  two  interior 
sources  and  for  an  externally  illuminated  living  mouse. 
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X  an  X.  ram  X,  mm 

a  b  c 

Fig-7  Surface  reconstruction  by  ML  method  for  phantom  with  two  interior  sources  (a),  3D- 
reconstruction  for  sources  in  air  with  separation  2.52  mm  (b)  and  3D  surface  (skin)  reconstruction 
of  nude  mouse  from  reflected  light  (c). 


The  volumetric  reconstruction  of  the  source  is  a  very  important  part  of  LET,  and  this  problem 
is  formulated  as  an  inverse  source  problem  based  on  diffusion  approximation42  44  .  The  solution 
uniqueness  for  BLI  reconstruction  can  be  assured  by  incorporation  of  a  priori  information  on  the 
animal  anatomy  and  optical  properties  24,45-47 .  Good  results  for  3D-  optical  reconstruction  have  also 
been  obtained  by  using  the  Finite  Element  Method48  52 .  In 50  two  sources  with  ~10  mm  separation 
were  reconstructed  using  an  FEM  algorithm.  We  will  show  below  our  reconstruction  results  for  the 
same  separation. 

We  present  another  technique  -  the  iterative  deblurring  EM  method53-55  based  on  diffusion 
approximation.  In  our  approach  good  reconstruction  can  be  obtained  if  the  sources  intensities  on  the 
surface  permit  some  separation  (see  Fig.6).  Our  experimental  data  in  this  case  for  backprojection 
onto  x  -axis  follows  this  (see  Fig.  8),  and  we  expect  a  good  reconstruction. 


Fig.8  Backprojection  histograms  for  1%  Intralipid  phantom  with  two  sources  located  in  0  mm 
and  10.16  mm  (a)  and  for  phantoms  with  different  source  positions  on  x-axis  in  inches  (b).  On  y- 
and  z-axis  sources  located  in  0  mm  position  for  each  phantom.  Voxel  size  =  0.5  mm. 
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We  can  take  into  account  reflection  on  boundary33-37  estimate  values  for  each  voxel,  and 
subtract  this  reflection  contribution.  Next  for  each  voxel  we  can  determine  a  “zero-th”  order 
approximation  (i.e.  initial  guess  as  a  input  image)  using  the  solution  for  point  source  photon  flux: 

<P°i  =  (P®<P  +  n)j  =  ~t~ - IWI  exP (-Mejrrjs)/ris  (13) 

Y.xVP'i  '=1 

i= I 

where  p  is  a  point-spread  function,  (p  is  a  real  value  of  photon  flux,  h  is  a  noise  term,  and 


xtj  is  the  path  length  of  line  /  through  voxel  j,  p\  and  are  intensities  in  CCDs  from  source 

in  media  and  in  air,  and  S  is  number  of  sources.  After  correcting  all  voxels  for  boundary  conditions 
we  apply  the  deblurring  EM  algorithm  for  3D-image  reconstruction: 


p® 


<Pj 


p®<pn. 


(14) 


where  ®  denotes  convolution  procedure  based  on  the  discrete  Fast  Fourier  Transform56 . 
Parameters  p  and  p  are  the  deblurring  kernels  (for  example  Gaussian  functions)  and  can  be 
changed  with  number  of  iterations  n  . 
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Key  Research  Accomplishments 

•  Gantry  and  camera  housings  have  been  designed,  built  and  tested  . 

•  Camera  and  support  motion  assemblies  were  built. 

•  The  CCDs  and  associated  first  generation  optical/scintillator  components  were  built  and 
tested. 

•  The  data  acquisition  system  has  been  completed  and  assembled  with  the  above  components. 
These  have  been  used  with  a  test  optical  imaging  system  based  on  SITe  CCD  chips,  now 
completed  and  in  use. 

•  The  5-head  system  was  constructed,  but  we  experienced  CCD  chip  failure.  As  explained 
below,  this  caused  us  to  move  to  a  new  CCD  assembly,  with  the  same  price  but  different 
characteristics.  This  required  a  redesign  of  the  electronics,  still  in  progress. 

•  Single  camera  tests  were  successful,  showing  high  sensitivity. 

•  Development  of  optical  MLEM  reconstruction  software  has  progressed,  and  we  have 
reported  our  first  results  in  tomographic  imaging  of  light  emission.  Our  images  show  that  it 
is  possible  to  detect  small  lesions  not  only  subcutaneously  but  also  at  depth  in  vivo,  and  we 
characterizef  the  advantages  of  this  method  for  quantitative  imaging.  3D  images  were 
obtained  and  reported  at  the  ERA  of  HOPE  symposium. 

•  Phantom  testing  has  been  completed  including  both  homogeneous,  cylindrical,  and 
heterogeneous,  mouse-shaped  phantoms.  We  have  also  demonstrated  micro-SPECT  (Single 
Photon  Emission  Computed  Tomography)  imaging  at  high  resolution. 


Reportable  Outcomes 

1)  ERA  of  HOPE  abstract  P64-3 

INVESTIGATION  OF  METASTATIC  BREAST  TUMOR  HETEROGENEITY  AND  PROGRESSION 
USING  DUAL  OPTICAL/SPECT  IMAGING 

Peter  P.  Antich,  Matthew  A.  Lewis,  Edmond  Richer,  Nikolai  Slavine,  Todd  Soesbe,  Xiufeng  Li, 
Anca  Constantinescu,  Allen  Harper,  Ralph  P.  Mason 

The  immediate  goal  of  our  project  is  to  image  tumor  growth,  metastatic  development 
and  vascular  changes,  both  to  characterize  tumor  dynamics  during  growth  for  application 
in  diagnostic  and  prognostic  imaging,  and  to  aid  in  the  development  of  new  therapies  by 
dissecting  the  effects  of  treatment.  More  specifically,  our  technical  goal  for  the  first  year  is 
to  develop  he  device  necessary  to  image  by  nuclear  medicine  and  optical  techniques.  We 
have  begun  assembling  the  physical  and  biological  components  of  our  research. 

Optical  imaging  techniques  have  been  developed  to  characterize  cells  and  molecules 
in  the  laboratory.  They  are  currently  used  with  microscope  techniques  in  excised  tissues  or 
living  cells.  They  are  now  being  extended  to  the  analysis  of  biologic  phenomena  in  living 
mammalian  organisms  and  have  found  a  strong  application  in  drug  screening.  Their  first 
applications  to  noninvasive  human  studies  are  envisaged  in  the  detection  of  oxygen 
saturation.  We  are  interested  in  advancing  this  technology  to  obtain,  e.g.,  oxygen- 
dependent  3D  images.  The  successful  translation  to  the  clinic  depends  upon  a  crucial 
question  -  can  small  lesions  be  detected  at  depth  or  is  this  method  limited  to  contact  or 
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small  distance  imaging?  We  had  previously  completed  a  three-dimensional  imager  of  light 
emitted  via  bioluminescence  both  in  its  hardware  and  image  reconstruction  software.  We 
report  here  on  initial  results  in  tomographic  imaging  of  light  emission.  We  present  images 
showing  that  it  is  possible  to  detect  small  lesions  not  only  subcutaneously  but  also  at  depth 
in  the  mouse  lungs  and  viscera,  and  characterize  the  advantages  of  this  method  for 
quantitative  imaging. 

We  have  also  demonstrated  micro-SPECT  imaging  of  bone,  in  order  to  assess  our 
ability  to  detect  vascular  phenomena  and  bone  metastases  at  high  resolution.  Nuclear 
medicine  is  a  well  accepted  clinical  tool,  for  which  agents  and  methods  exist  reportedly 
able  to  detect  and  potentially  stage  disease  in  the  clinical  environment.  We  report  here  on 
our  initial  progress  in  SPECT,  and  on  the  advantages  of  tomography  initially  for  small 
animals,  showing  submillimeter  resolution  in  a  mouse  when  imaging  bone.  This  form  of 
imaging  may  be  used  to  detect  bone  metastases.  POSTER  ATTACHED 

2J 

NOVEL  APPLICATION  OF  BIOLUMINESCENT  IMAGING:  INTERROGATING  ACUTE 
EFFECTS  OF  THE  VASCULAR  TARGETING  AGENT  COMBRETASTATIN. 

Dawen  Zhao,  Allen  Harper,  Edmond  Richer,  Li  Liu,  Nikolai  Slavine,  Jerry  W.  Shay,  Peter 
P.  Antich  and  Ralph  P.  Mason 

Cancer  Imaging  Program,  Department  of  Radiology  and  Cell  Biology,  UT  Southwestern, 
Dallas,  TX,  USA  75390 

Vascular  targeting  agents  promise  therapeutic  efficacy  by  selectively  cutting  of  the  blood 
supply  to  tumors  inducing  massive  cell  death.  Combretastatin  (CA4P)  represents  a  lead 
drug  and  is  in  advanced  clinical  trials.  Numerous  reports  indicate  acute  changes  in 
perfusion  observed  by  MRI  reporting  rapid  vascular  shutdown  in  tumors.  We  recently 
measured  acute  hypoxiation  within  90  mins  following  CA4P  administration  to  rats  with 

breast  tumors  using  I^F  MRI  oximetry  (Zhao  et  al.,  IJROBP  2005).  We  have  now  applied 
BLI  to  probe  the  acute  effects  of  CA4P  on  human  breast  231  adenocarcinomas  transected 
to  stably  express  luciferase. 

231  cells  were  transfected  with  lentivirus  to  express  luciferase  and  highly  expressing 

clones  isolated.  10®  cells  were  implanted  in  the  flank  of  nude  mice  and  allowed  to  grow  to 
~  6  mm  diameter.  Sty  this  stage  mice  were  anesthetized  (isoflurane/C>2).  A  solution  of 

luciferin-D  (sodium  salt  450  mg/kg)  was  administered  SC  in  the  neck  region  and  light 
images  acquired  immediately  using  one  camera  of  our  Light  Emission  Tomography 
System  (LETS).  Serial  images  (30  s  each)  were  acquired  over  a  period  of  20  -  30  mins 
and  the  light  intensity  -  time  curves  evaluated.  Saline  or  CA4P  in  saline  (120  mg/kg)  were 
injected  IP  immediately  after  baseline  BLI  and  then  2  h  and  24  h  later  the  BLI  time  course 
was  repeated.  We  also  undertook  3D  imaging  by  acquiring  multiple  images  simultaneously 
using  3  cameras  arranged  in  a  circular  gantry  around  the  mouse.  Data  were  reconstructed 
using  a  statistical  MLEM  approach  and  deconvolution  to  provide  fully  3D  images. 

Control  tumors  showed  intense  light  emission  peaking  within  8  mins  and  generally 
decreasing  to  about  50-70%  after  20  mins.  Repeat  measurements  after  saline  infusion 
essentially  indicated  similar  kinetics.  By  contrast,  CA4P  led  to  a  significantly  lower  light 
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emission  (peak  ~  2  to  10  times  lower)  and  delayed  peak  emission.  24  h  later  signal 
remained  considerably  decreased.  Traditional  BLI  provides  a  planar  image  only,  but 
tumors  are  3D  and  known  to  exhibit  heterogeneity  particularly  after  vascular  targeting 
agents.  LETS  was  successful  providing  3D  representation  of  the  tumors. 

While  BLI  has  been  reported  extensively  in  terms  of  measuring  tumor  burden  with  respect 
to  growth  and  response  to  therapy,  we  believe  that  application  to  examine  acute  effects  of 
vascular  targeting  agents  is  a  novel  application.  It  can  provide  a  cheap  high  throughput 
method  for  evaluating  novel  drugs  and  drug  combinations  and  scheduling. 

Supported  by  the  NCI  Cancer  Imaging  Program  pre-ICMIC  CA86354,  DOD  Breast 
Cancer  Initiative  (DAMD  170310363)  (DZ)  and  Simmons  Cancer  Center. 


3)  University  of  Texas  Molecular  Medicine  Symposium  Invited  Speaker 

Three-Dimensional  Imaging  of  Cellular  and  Molecular  Events  in  the  Living 
Animal  by  Light  Emission  Tomography 

The  emerging  field  of  in  vivo  “molecular  imaging”  encompasses  a  diverse  array  of 
noninvasive  imaging  methods  which  enable  gene  expression  and  other  molecular  biology 
assays  to  be  detected  and  visualized  in  living  animals  Cell  assays  have  been  developed 
to  study  genetic  regulatory  elements,  protein  trafficking  and  protein  function.  Extension  of 
this  type  of  analysis  to  intact  living  animals,  and  eventually  translation  to  humans,  will 
enhance  our  understanding  of  gene  regulation  and  protein  function  in  functional  tissues 
and  organs.  Molecular  imaging  is  now  performed  with  MRI,  PET,  SPECT  and  CT,  as  well 
as  by  the  latest  light  emission  imaging  techniques,  in  particular  bioluminescence  and 
fluorescence  imaging. 

When  sets  of  images  obtained  by  different  techniques,  or  different  reporter  molecules, 
measure  and  display  different  properties,  it  is  of  great  interest  to  analyze  their  interrelation. 
To  do  so,  it  is  preferable  to  achieve  spatio-temporal  coregistration  of  data  sets  obtained 
with  the  different  methods.  To  this  end,  molecular  imaging  methods  should  preferably  be 
three-dimensional  and  also  quantitative,  so  that  the  data  extracted  from  its  images  can  be 
used  to  reveal  mechanisms  through  which  the  biological  processes  unfold. 

This  talk  focuses  on  techniques  and  methods  to  perform  quantitative,  three- 
dimensional  bioluminescent  imaging  (BLI).  The  technique  is  based  on  transfection  of  the 
luciferase  gene  (North  American  firefly,  Photinus  pyralis),  which  is  readily  introduced  into 
diverse  cell  lines  with  selection  for  most  effective  expression.  When  D-luciferin  is 
administered  to  the  cells  or  to  an  animal  host,  a  bioluminescent  reaction  emits  light  that 
can  be  detected.  Our  current  research  focus  is  on  imaging  the  evolution  and  treatment  of 
primary  tumors  and  distant  metastases.  Detection  of  metastatic  foci,  study  of  the  dynamics 
of  metastatic  development,  and  study  of  response  to  therapy  of  primary  tumor  and 
metastatic  foci  are  tasks  which  demand  high  sensitivity  and  resolution  in  the  context  of 
noninvasive  whole-body  imaging,. 

To  date,  BLI  has  been  mostly  limited  to  a  planar  technique,  where  the  effectiveness  of 
the  2D  method  is  restricted  in  terms  of  localization  of  cell  population  (visualization  at  depth) 
and  size  estimation  (quantitation  of  tumor  burden  and  progression).  Due  to  light  scattering 
and  absorption,  3-D  reconstruction  of  light  emission  from  a  source  at  depth  poses  major 
difficulties.  The  intensity  and  spatial  distribution  of  the  BLI  signal  are  strongly  influenced  by 
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the  intervening  tissues,  which  can  be  accounted  for  numerically  if  the  anatomic  distribution 
of  tissues  and  their  optical  properties  are  known. 

Here,  we  describe  an  imaging  method  to  accomplish  this  task,  light  emission 
tomography  (LET),  obtained  by  applying  to  data  obtained  all  around  the  mouse  an 
algorithm  based  on  image  reconstruction  techniques  well-established  for  single-photon 
emission  tomography,  but  differing  in  dynamic  factors  which  account  for  photon  absorption 
and  diffusion.  The  reconstruction  method  is  validated  in  simple  phantoms  of  appropriate 
dimensions,  showing  its  ability  to  overcome  diffusion  effects  and  to  achieve  optimal 
quantitative  power.  Representative  images  obtained  in  living  animals  are  shown  to 
document  the  ability  of  the  method  to  image  small  lesions  at  depth. 

In  summary,  three-dimensional  (3D)  optical  imaging  of  bioluminescence  has  the  potential 
to  provide  unique  insight  into  tumor  development  and  metastasis.  Compared  to  planar  BLI, 
it  enhances  sensitivity  and  detectability.  The  method  permits  a  quantitative  assessment  of 
cancer  cell  burden  in  living  animals,  providing  information  both  on  cell  numbers  and, 
combined  with  anatomical  information,  their  localization  in  organs. 

4) 

Richer  E,  Slavine  N,  Lewis  MA,  Tsyganov  E,  Gellert  GC,  Dikmen  ZG,  Bhagwandin  V,  Shay  JW, 
Mason  RP,  Antich  PP:  Three  Dimensional  Light  Emission  Tomography  using  Multiple  Rotating 
CCD  Cameras.  Third  Annual  Meeting  of  The  Society  for  Molecular  Imaging,  9-12  September, 

2004,  St.  Louis  POSTER  ATTACHED 

5) 

Richer  E,  Lewis  MA,  Smith  B,  Li  X,  Seliounine  S,  Mason  RP,  Antich  PP:  Comparison  of  CsI(Tl) 
and  Scintillating  Plastic  in  a  Multi-Pinhole/CCD-based  Gamma  Camera  for  Small  Animal  Low 
Energy  SPECT.  Small-animal  SPECT  Imaging:  Proceedings  of  The  University  of  Arizona  Center 
for  Gamma-Ray  Imaging  Workshop  on  Small-Animal  SPECT,  Chapter  13,  Springer-Verlag,  2005. 

PREPRINT  ATTACHED 

6) 

Lewis  MA,  Arbique  G,  Richer  E,  Slavine  N,  Jennewein  M,  Constantinescu  A,  Brekken  R,  Guild  J, 
Tsyganov  EN,  Mason  RP,  Antich  PP:  Projection  and  pinhole  based  data  acquisition  for  small 
animal  SPECT  using  storage  phosphor  technology.  Small-animal  SPECT  Imaging:  Proceedings  of 
The  University  of  Arizona  Center  for  Gamma-Ray  Imaging  Workshop  on  Small-Animal  SPECT, 
Chapter  26,  Springer-Verlag,  2005.  PREPRINT  ATTACHED 

7; 

Z.G.  Dikmen,  G.  Gellert,  P.  Dogan,  R.  Mason,  P.  Antich,  E.  Richer,  W.E.  Wright,  J.W.  Shay,  “ 

A  New  Diagnostic  System  in  Cancer  Research:  Bioluminescent  Imaging",  Turk.  J.  Med.  Sci.  35,  65- 
70  (2005). 

5; 

Paroo  Z,  Bollinger  RA,  Braasch  DA,  Richer  E,  Corey  DR,  Antich  PP,  Mason  RP  (2004) 
Validating  bioluminescence  imaging  as  a  high-throughput,  quantitative  modality  for 
assesing  tumor  burden.  Molecular  Imaging  3(2):  1 17-124. 
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Conclusions 


Successful  completion  of  the  instrument  development  specific  aims  in  Year  1  will  allow  the 
initiation  of  experiments  in  the  animal  model  that  were  proposed  in  Year  2.  Initial  imaging 
experiments  in  tumor-bearing  mice  are  in  progress  at  present. 
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1.  INTRODUCTION 

SPECT  with  submillimeter  resolution  (/uSPECT)  is  needed  for  small 
animal  imaging  and  differs  from  clinical  SPECT  -  a  rich  and  well  es¬ 
tablished  methodology  of  nuclear  medicine  -  in  that  its  focus  is  on  the 
specific  requirements  of  high  resolution  imaging  over  small  distances 
with  matching  sensitivity,  requiring  differences  in  collimator  schemes 
and  image  reconstruction  methods  specialized  for  that  task.  Current 
efforts  towards  this  goal  have  focused  on  adapting  established  technolo¬ 
gies  from  clinical  to  animal  imaging.  Currently  available  devices  have 
resolution  in  the  range  of  1-3  mm  and  adequate  sensitivity,  but  we  be¬ 
lieve  that  these  parameters  should  be  improved  to  facilitate  molecular 
imaging.  Prototypes  are  being  developed  [Ogawa  et  al.,  1997;  Ogawa. 
et  al.,  1998;  MacDonald  et  al.,  2001;  McElroy  et  al.,  2001;  Schramm 
et  al..  2001]  in  response  to  a  rising  level  of  interest  by  pharmaceutical 
companies  and  molecular  biologists  (Russo,  1998). 

Nevertheless,  a  significant  factor  for  penetration  of  /zSPECT  into  gen¬ 
eral  usage  remains  the  cost  and  ease  of  use  of  a  system  with  adequate 
resolution  for  studies  in  mice  and  rats.  In  the  case  of  small  animal  imag- 
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ing,  pinhole  collimator  and  a  closed  detector-object,  configuration  can 
be  used  with  SPECT.  This  presents  the  advantage  to  minimize  the  time 
for  data  collection,  allowing  a  superior  detector  sensitivity  at  high  spa¬ 
tial  resolution  by  using  multiple  pinholes  and  crystals.  A  pinhole  size  of 
0.5  mm  would  provide  the  sub-millimeter  resolution  required  to  study 
small  animals  such  as  mice. 

2.  CCD-BASED  GAMMA  CAMERA 

With  the  goal  of  achieving  high  resolution  with  appropriate  sensitiv¬ 
ity  for  use  in  small  animal  studies,  we  have  developed  a  multi-pinhole 
gamma  camera  based  on  thallium-doped  cesium  iodide  [CsI(Tl)]  crystals 
and  wavelength-matched  charge-coupled  devices  (CCDs). 

A  Scientific  Imaging  Technologies  Inc.  (SITe)  SIA502AB1  CCD,  a 
512x512  pixels,  non-color,  back-illuminated,  full  frame  image  sensor  was 
selected  due  to  its  excellent  sensitivity  (greater  than  85  %  quantum 
efficiency  from  400  nm  to  750  nm).  The  pixel  size  is  24  x  24  fim  providing 
a  large  well  capacity  of  350,000  e".  Coupled  with  a  good  sensitivity  (2.6 
low  dark  current  (20  pA/cm2  at  20°C),  and  low  readout  noise 
(5  e“)  tills  provides  a  dynamic  range  of  70,000,  more  than  full  16  bits. 

The  CCD  is  be  incorporated  into  a  self-contained,  cooled  camera, 
equipped  with  electronic  circuitry  and  optics.  The  principal  method  for 
noise  suppression  with  CCD  technology  is  cooling  the  chip  itself.  In 
our  proposed  implementation,  the  CCD  camera  employs  a  five  stage 
thermoelectric  cooling  element  equipped  with  a  controller  circuit  that 
keeps  the  operational  temperature  of  the  CCD  at  -40°C.  An  external 
heat  exchange  unit  eliminates  the  generated  heat.  The  net  effect  will  be 
a  substantial  reduction  of  the  dark  noise  generated  inside  the  CCD,  and 
thus  better  image  quality  and  longer  possible  exposures.  The  internal 
chamber  of  the  camera  is  kept  under  vacuum  to  allow  better  cooling  and 
eliminate  ice  crystals  formation  on  the  sensitive  element. 

The  light  exposure  is  controlled  using  a  mechanical  shutter.  The 
shutter  allows  for  short  exposure  times,  down  to  5  msec,  which  are  not 
attainable  when  an  electronic  shutter  is  employed.  Also,  the  mechanical 
shutter  prevents  light  from  affecting  the  CCD  during  the  readout  phase, 
thus  eliminating  the  smear  effect  observed  in  full  frame  CCD  cameras 
without  a  shutter.  To  widen  the  dynamic  range,  the  CCD  readout  em¬ 
ploys  a  16-bit,  A/D  conversion  chip  equipped  with  a  low  noise  preampli¬ 
fier  ana  correlated  double  sampling  circuitry.  This  architecture  provides 
high  dynamic  range  combined  with  very  low  readout  noise.  Finally,  a 
flexible  pixel  binning  scheme  is  implemented  allowing  for  better  signal 
to  noise  ratio  in  the  extremely  low  light  level  conditions  that  expected 
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with  higher  resolution  imaging.  Optical  coupling  of  the  CCD  to  the 
scintillating  crystal  is  implemented  using  an  off-the-shelf  low  reflectance 
optical  lens. 

We  selected  CsI(Tl)  as  the  scintillator  rather  than  the  more  traditional 
gamma  camera  scintillator,  Nal(Tl).  Nal(Tl)  has  a  scintillation  wave¬ 
length  of  maximum  emission  at  415  nm,  ideally  matched  to  the  efficiency 
profile  of  photomultiplier  tubes  (PMTs),  which  have  a  typical  quantum 
efficiencies  of  20%.  CsI(Tl),  on  the  other  hand,  has  a  maximum  emission 
at  550  nm,  which  is  well  into  the  region  where  solid-state  detectors  are 
dominant.  While  most  scintillator  property  tables  list  CsI(Tl)  as  having 
45%  the  photoelectron  yield  as  compared  to  Nal(Tl),  this  comparison  is 
deceptive  as  it  holds  true  only  for  PMTs;  using  CCDs,  the  light  yield  of 
Csl  is  53%  higher  than  Nal.  Further  advantages  are  sturdier  mechanical 
properties,  limited  hygroscopicity,  and  higher  density.  This  last  factor 
will  allow  thinner  scintillators  with  the  same  capture  efficiency. 

Two  potential  drawbacks  exist  for  this  gamma  camera  design,  intrinsic 
to  the  scintillator  and  the  detector.  First,  CsI(Tl)  is  a  slow  scintillator 
(900  ns  versus  230  ns  for  Nal(Tl)).  However,  this  is  not  a  problem  since 
timing  and  triggering  is  not  needed  because  the  CCD  is  an  integrating 
device,  and  subsequently  an  exposure  will  contain  the  recorded  detection 
of  many  high  energy  photons.  Second,  and  more  importantly,  energy 
discrimination  is  not  possible.  Two  facts  make  its  application  to  mouse 
imaging  possible  and  free  of  excessive  scatter  contributions.  First,  our 
major  applications  of  the  proposed  SPECT  is  to  1-125  imaging,  in  which 
case  the  existence  of  a  K-edge  at  32.2  keV  for  the  crystal  is  sufficient 
to  accept  the  35  keV  gamma  emitted  by  this  nuclide  and  reject  lower 
energy  noise.  Second,  the  completed  system  is  not  intended  for  use  in 
humans,  where  gamma  rays  have  an  appreciable  probability  of  scatter 
before  exiting  the  body,  but  for  a  small  animal,  so  that  scatter  rejection 
is  less  important  and  can  be  accounted  for  with  absorption  corrections. 

The  intrinsic  resolution  of  the  gamma  detector  was  determined  using 
a  100  fm i  slit  collimator  displaced  across  the  surface  of  the  scintillating 
crystal  using  a  micrometric  positioning  device.  Figure  (13.1)  shows  the 
image  from  a  1-125  line  source  displaced  0.5  mm  across  crystal  face.  The 
corresponding  profile  is  presented  in  Fig.  (13.2).  The  intrinsic  resolution 
of  the  camera  was  determined  to  be  below  0.4  mm  (see  Antich  et  al., 
2003).  Considering  a  field  a  view  of  50  mm,  suitable  for  an  animal  the 
size  of  a  mouse,  the  number  of  pixels  required  to  record  an  image  in  a 
non  magnifying  geometry  is  less  than  128.  Compaxing  with  the  much 
higher  resolution  of  the  CCD  (512x512  pixels)  one  concludes  that  a  set 
of  up  to  16  pinhole  projections  can  be  acquired  on  a  single  CCD  image. 
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Figure  13.1  100  fim  1-125  line  source 

displaced  0.5  mm  across  crystal  face. 


! 


Pixel* 


Figure  13.2.  Profile  of  the  line  source 
displaced  0.5  mm  across  crystal  face. 


Thus,  in  out  design,  collimation  is  provided  by  a  3x3  array  of  pinholes 
that  provide  some  field-of-view  overlap  in  the  transverse  plane,  but  are 
non- degenerate  in  the  longitudinal  axis  of  the  animal.  A  schematic  ren¬ 
dering  of  the  geometry  of  this  set-up  and  an  actual  multiple-pinhole 
image  is  shown  in  Figs.  (13.3)  and  (13.4). 


Figure  13.3.  Field  of  view  of  pinhole  Figure  IS.4.  Convergent  9-piuhole 

collimators.  (1.5  mm)  collimator,  2  400  /iCi 

(14.8  MBq)  line  sources,  30  minute 
acquisition. 
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3.  PLASTIC  SCINTILLATORS 

For  imaging  with  low  energy  iodine-125  labeled  agents  (<35.5  keV), 
we  further  consider  scintillating  plastic  as  an  alternative  to  CsI(Tl). 
Plastic  scintillators  have  the  advantage  of  being  mechanically  sturdier 
than  crystals,  have  less  hygroscopicity,  and  more  importantly  can  be 
shaped  in  any  desired  shape  and  dimension.  Moreover,  the  light  emis¬ 
sion  can  be  tuned  to  the  optimum  wavelength  sensitivity  of  the  CCDs 
using  waveshifting  doping  materials.  A  potential  drawback  is  their  lower 
density,  which  require  thicker  scintillators  to  achieve  a  good  detector 
cross-section  with  negative  effects  on  resolution. 

In  order  to  determine  the  efficiency  of  the  plastic  scintillators  we  com¬ 
pared  the  light  output  produced  by  a  125-1  line  source  (400  pCi)  in  direct- 
contact  with  a  5.6  mm  BC-408  scintillator  (Bicron  Inc.)  and  a  2  mm 
Cil(Tl)  crystal.  The  light  output  normalized  to  the  thickness  shows  a 
factor  of  237  in  favor  of  the  CsI(Tl)  crystal. 

A  set  of  PVT  based  scintillators  containing  different  microcrystals  in 
two  concentrations  were  also  tested  with  the  results  presented  in  Table 
(13.1).  The  concentration  is  given  in  mg  of  microcrystal  material  per  g 
of  PVT,  and  the  light  output  in  arbitrary  light  units  per  minute  of  ex¬ 
posure.  (ALU /min)  as  obtained  from  the  CCD  camera.  The  light  output 
normalized  to  the  scintillator  material  thickness  is  also  presented.  All 
tested  microcrystals  significantly  improved  the  efficiency  of  the  scintil¬ 
lator,  with  LSO  in  100  mg/g  concentration  showing  spectacular  results, 
increasing  the  thickness  normalized  light  output  almost  600  times  com¬ 
pared  to  the  PVT  alone.  Comparing  the  PVT  containing  LSO  with  the 
CsI(Tl)  crystal  we  observed  that  while  the  total  light  output  is  larger 
for  the  crystal  by  a  factor  of  2.26  the  normalized  light  output  is  higher 
for  the  PVT  containing  LSO  by  a  factor  of  1.58. 


Table  13.1.  Comparison  of  plastic  scintillators  light  output  with  different  microcrys- 
tais. 


Microcrystal 

Material 

Microcrystal 

Concentration 

Thickness 

[mmj 

Light  Output 
[ALU /min] 

Normalized 
Light  Output 

ErSO 

50 

0.5 

375 

750 

ErSO 

100 

0.5 

625 

1,250 

YbSO 

50 

0.35 

525 

1,500 

YbSO 

100 

0.35 

888 

2,537 

LSO 

50 

0.35 

10,125 

28,929 

LSO 

100 

0.35 

16,065 

45,900 

None 

0 

5.6 

431 

77 
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The  intriasic  resolution  of  the  scintillating  plastic  LSO/CCD-coupled 
gamma  camera  was  determined  by  the  same  method  to  be  <0.5  mm, 
comparable  to  our  previously  reported  results  for  the  CsI(Tl)  based  de¬ 
tector.  An  image  of  2  400  /iCi  (14.8  MBq)  line  sources,  obtained  using 
LSO  100  mg/g  scintillator  and  a  1.5  mm  pinhole  collimator  is  presented 
in  Fig.  (13.5). 


Figzm  13.5.  Pinhole  (1.5  mm)  collimator.  2  400  /*=Ci  (14.8  MBq)  line  sources,  30 
minute  acquisition. 


A  potential  problem  with  the  natural  radioactive  background  of  176Lu 
in  LSO  has  been  identified  with  single  photon  imaging  [Huber  et  al., 
2002].  However,  since  the  natural  background  radiation  spectrum  is 
primarily  above  300  keV,  the  thin  section  of  the  doped  plastic  scintillator 
will  minimize  unwanted  noise  from  this  source. 

4.  CONCLUSIONS  AND  FUTURE  WORK 

We  have  made  considerable  progress  in  designing  and  constructing 
devices  which  will  allow  us  to  compare  quantitatively  different  designs 
and  technologies  for  high  resolution,  high  sensitivity  ^.SPECT.  In  future 
studies  we  will  optimize  with  respect  to  the  resolution  and  sensitivity 
the  plastic  scintillator  thickness  and  the  microcrystal  concentration. 

Challenges  also  remain  in  the  development  of  the  image  reconstruction 
algorithms  suitable  for  multiple-pinhole  CCD-based  /zSPECT  cameras. 
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Chapter  26 

Projection  and  Pinhole  based  Data  Acquisition 
for  Small  Animal  SPECT  using  Storage  Phosphor 
Technology 


M  A  Lewis,  G  Arbique,  E  Richer,  N  Slavine,  M  Jermewein*,  A  Con- 
stantinescu,  R  Brekken,  J  Guild,  E  N  Tsyganov,  R  P  Mason,  and  P  P 
Antich 

The  University  of  Texas  Southwestern  Medical  Center  at  Dallas 

1.  Introduction 

Three-dimensional  Single  Photon  Emission  Computed  Tomography 
(SPECT)  can  provide  high-resolution  insight  into  biomolecular  distri¬ 
bution  and  pharmacokinetics.  However,  instrument  availability  and 
distribution  is  limited  at  present,  and  imaging  times  can  be.  consider¬ 
able.  To  evaluate  the  large  array  of  novel  agents  which  are  becoming 
available,  we  find  that  storage  phosphor-based  in  vivo  imaging  can  pro¬ 
vide  an  important,  rapid-throughput  transition  from  the  traditional  ex 
vivo  sacrifice/gamma-counting  and  autoradiography  to  full  time-course 
SPECT. 

Storage  phosphor  (SP)  technology  has  found  wide-spread  usage  in 
both  the  radiology  clinic  and  in  the  molecular  biology  laboratory.  With 
a  linear  sensitivity  latitude  exceeding  that  of  plain  film,  photostimula- 
ble  luminescence  (PSL)  screens  can  facilitate  high-throughput  screening 
and  autoradiography  for  a  variety  of  small  animal  models.  Despite  the 
impressive  growth  and  improvement  in  small  animal  nuclear  imaging, 
the  determination  of  appropriate  time  and  duration  of  imaging  for  a 
given  animal  model  and  imaging  agent  is  still  crucial  in  study  design 
and  analysis.  Storage  phosphor  technology  is  well-suited  for  addressing 
this  issue,  either  by  producing  quantitative  planar  projection  or  low- 
resolution  tomographic  images.  With  7-ray  stopping  cross-sections  that 
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are  not  radically  different  from  film/screen  systems,  the  phosphor  image 
plate1  is  well-suited  for  low  energy  gamma  detection,  especially  for  imag¬ 
ing  studies  with  iodine-125  (125I)  where  tissue  scattering  is  minimized. 

2.  Background 

2.1  Storage  phosphor  technology 

Photostimulable  storage  phosphor  (PSP2)  screens  were  developed  and 
introduced  in  the  1980s  [Sonoda  et  al.,  1983]  as  a  replacement  for  film/screen 
x-ray  detection  systems.  Today,  SP  technology  is  the  key  component  in 
commercial  film-less  imaging  systems  grouped  under  the  generic  term 
Computed  Radiography  (CR)  [Arbique  et  al.,  2003].  This  technology 
is  distinct  from  so-called  Digital  Radiography,  where  x-ray  detection 
is  digital  from  the  moment  of  detection.  Subsequently,  SP  technology 
was  adopted  by  most  suppliers  of  autoradiography  equipment  [Johnston 
et  al.,  1990].  While  the  average  molecular  biologist  today  does  not  have 
ready  access  to  a  small  animal  SPECT  system,  he  or  she  most  likely 
possesses  or  has  available  an  SP  scanner  or  reader  that  is  routinely  used 
for  gel  and  blot  applications.  Compared  to  the  standard  CR  reader  in  a 
radiology  clinic,  these  scanners  typically  have  higher  spatial  resolutions 
(down  to  50  /im)  and  greater  digital  resolution  (16  vs.  12  bit).  Although 
the  current  trend  is  towards  more  expensive  scanners  that  can  also  read 
direct-  or  chemi-fluorescence  assays,  a  multipurpose  SP  system  can  be 
obtained  for  a  fraction  of  the  cost  of  a  minimal  small  animal  SPECT 
system. 

A  full  physical  description  of  SP  technology,  including  image  per¬ 
formance  parameters  (MTF,  NPS,DQE)  can  be  found  elsewhere  [von 
Seggern,  1999;  Rowlands,  2002].  In  summary,  the  interaction  of  a  high- 
energy  photon  with  a  BaF(Br,I):Eu2+  phosphor  grain  generates  electrons 
and  holes  that  are  trapped  in  the  immediate  fluorohalide  matrix  area. 

A  latent  image  is  formed  in  proportion  to  the  density  of  the  trapped 
carriers.  Upon  subsequent  photostimulation  using  a  red  light  source 
(typically  a  laser  at  590-680  nm),  the  electron  is  liberated  from  its  trap 
and  recombines  with  a  hole.  The  energy  generated  in  this  process  is 
resonantly  transfered  to  a  doped  Europium  ion,  which  then  decays  with 
characteristic  luminescence  (390  nm).  This  luminescence  is  captured  us¬ 
ing  a  photomultiplier  tube  and  digitized.  With  an  appropriate  reader,  a 
typical  SP  image  plate  will  have  a  linear  dynamic  range  up  to  a  thousand 
times  greater  than  film,  spanning  4-5  orders  of  magnitude,  and  reusable 
for  tens  of  thousands  of  exposures.  Due  to  thermal  stimulation,  the 
latent  image  does  fade  (decay  constant  approximately  9  hours),  so  im¬ 
mediate  reading  of  PSL  screens  is  advisable. 
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2.2  Low  energy-integrating  SPECT 

The  hypothesis  of  low  energy-integrating  SPECT  is  that  reasonable 
images  can  be  obtained  using  7-ray  detectors  that  integrate  activity 
counts  and  have  no  energy  discrimination  for  scatter  rejection.  Monte 
Carlo  simulations  [Tenney,  2002]  have  shown  that  for  125I  in  a  mouse¬ 
sized  cylinder  of  water  the  scatter  fraction  remains  below  40%  and 
the  detective  quantum  efficiency  (DQE)  is  in  the  reasonable  range  of 
0.6  —  0.7,  which  can  exceed  the  DQE  of  the  SP  screen  [Rowlands,  2002]. 
Although  dual  isotope  imaging  is  precluded,  energy-integration  imaging 
shifts  the  burden  of  work  to  image  reconstruction,  where  scatter  correc¬ 
tion  represents  a  significant  challenge,  but  is  an  active  area  of  research 
for  all  nuclear  medicine  modalities. 

A  CCD-based  gamma  camera  has  been  recently  proposed  as  an  effi¬ 
cient  energy-integrating  detector  for  small  animal  SPECT.  Although  a 
CCD-based  camera  will  in  general  be  less  expensive  than  the  traditional 
position-sensitive  PMT,  the  intrinsic  resolution  of  the  CCD  design  will 
be  significantly  better,  with  characteristic  resolutions  below  400  [An- 
tich  et  al..  2003].  At  these  low  intrinsic  resolutions,  parallel-hole  collima¬ 
tors  are  not  well-matched.  In  addition,  compared  to  a  non-magnifying 
pinhole  geometry,  the  packing  of  multiple,  static  projection  is  limited. 

For  lower  resolution,  high-throughput  screening,  we  find  that  SP  im¬ 
age  plates  and  screens  are  an  acceptable  alternative  to  CCDs  in  an 
energy-integrating  detector.  If  an  infrastructure  for  reading  SP  tech¬ 
nology  is  already  in  place,  then  the  relative  cost  for  detectors  is  further 
reduced,  with  no  sacrifice  in  intrinsic  resolution.  With  this  savings,  we 
believe  that  small  animal  SPECT  using  SP  technology  may  be  possible 
for  routine  low-resolution,  high-throughput  screening. 

3.  Prototype 

Recently,  a  prototype  system  using  CR  image  plates  for  co-registered 
1-125  planar  scintigraphy  and  x-ray  radiography  was  introduced  [Boone 
et  al.,  2003].  In  this  scheme,  half  of  the  image  plate  is  used  for  a  planar 
emission  image.  After  exposure,  the  image  plate  is  precisely  translated 
so  that  a  transmission  radiography  can  be  taken  with  the  other  half  of 
the  plate.  These  two  images  are  subsequently  co-registered  in  software, 
using  the  known  displacement. 

Our  approach  differs  in  that  we  have  utilized  custom-sized  SP  image 
plates  as  energy-integrating  detectors  in  a  low-resolution  small-animal 
pinhole  emission  computed  tomography  geometry.  As  shown  in  Figure 
26.1,  a  30  mm  inner  diameter/ 5  mm  thick  brass  cylinder  serves  as  a 
support  for  eight  1.5  mm  diameter  pinholes  (also  in  brass,  with  90° 
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Figure  26.1.  Prototype  murine  pinhole  emission  computed  tomography  system  using 
using  energy-integrating  storage  phosphor  technology 


cone  of  acceptance) ,  spaced  in  45°  increments  around  a  ring.  In  this 
geometry,  the  field-of-view  in  a  plane  intersecting  the  axis  of  the  cylinder 
is  equivalent  to  the  30  mm  diameter  of  the  cylinder.  A  strip  of  8  small 
SP  image  plates  is  attached  to  the  support  structure  at  fixed  points  of 
known  distance  from  the  pinhole  centers.  Each  detector  is  placed  10 
mm  from  the  corresponding  pinhole  center.  In  this  geometry,  reading  of 
the  SP  screen  at  100  /i m  resolution  will  produce  128x128  and  256x128 
projections  for  60°  and  90°  cone-of- acceptance  pinholes,  respectively. 
For  125I  energies  (<  35.5  keV),  brass  is  sufficient  for  collimation  (HVL 
^  80/xm). 

Eight  24  mm  by  50  mm  SP  screens  were  cut  from  a  25.2x30.3  cm 
Fuji  CR  ST-VN  Imaging  Plate  (Fuji  Photo  Film,  Tokyo).  Attachment 
holes  of  1.6  mm  diameter  were  milled  at  the  end  of  each  screen  for 
attachment  to  the  imaging  system.  The  imaging  plates  were  connected 
in  a  continuous  strip  using  black  insulating  tape. 

The  system  is  placed  in  a  light-tight  enclosure  with  the  phantom  or 
animal  placed  inside  the  cylinder.  After  exposure3,  the  SP  detector  is 
carefully  removed  and  transported  in  a  light-tight  envelope  to  the  SP 
reader.  A  SP  scanner  with  the  ability  to  read  arbitrary  size  screens, 
such  at  the  Molecular  Dynamics  Storm  (Amersham  Biosciences,  Piscat- 
away,  NJ),  is  preferred,  but  in  practice  other  systems  may  work.  Special 
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Figure  26.2.  Eight  pinhole  projections  from  two  1.85  MBq  12yI  capillaries  (0.8  mm 
inner  diameter).  15  minutes  integration  exposure.  Signal- to-noise  ratio  ^  10 
.  Heterogeneities  due  to  air  bubbles  in  capillaries  and  infiltration  of 
sealant. 


attention  to  the  relationship  between  the  detector  coordinate  system  and 
the  angular  position  of  the  detector  should  be  maintained,  as  different 
SP  readers  may  report  data  in  mirror  orientations. 

4.  Results 

To  evaluate  the  prototype  system,  two  non-uniform  1.85  MBq  125I 
capillaries  with  inner  diameter  of  0.8  mm  were  prepared.  The  line  sources 
were  placed  off-center  in  the  field-of-view  for  the  ring  of  pinholes,  and 
were  non-parallel  with  an  average  separation  of  8  mm.  The  sources  were 
imaged  for  15  minutes,  followed  by  prompt  reading  of  the  SP  image 
plates.  Eight  pinhole  projections  of  the  line  source  can  be  found  in 
Figure  26.2. 

Using  unexposed  regions  of  the  SP  image  plates,  background  noise 
levels  were  determined  for  windowing.  Since  the  attachment  holes  ap¬ 
pear  as  voids  on  the  raw  projection  images,  detector  coordinate  origins 
were  calculated  as  the  midpoint,  and  projections  were  corrected  for  ro¬ 
tation.  A  3D  image  (48x48x60,  0.5  mm  voxels)  was  reconstructed  from 
these  eight  projections  using  the  maximum-likelihood  estimation  maxi¬ 
mization  (MLEM)  algorithm  with  convolution  for  resolution  modelling 
[Reader  et  al.,  2002;  Zinchenko  et  al,  2003].  Figure  26.3  provides  a 
maximum  intensity  projection  (MIP)  of  the  reconstructed  volume  that 
is  consistent  with  the  projections  in  Figure  26.2.  For  a  non-magnifying 
pinhole  geometry,  the  standard  pinhole  collimator  resolution  expression 


Rc  = 


a  +  6 

Ug 


a 


(26.1) 
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Figure  26.3.  Maximum  intensity  pro¬ 
jection  for  reconstructed,  heterogenous 
line  sources 


Figure  26.4 .  Cross  section  through 
line  sources  with  expected  broadening 
due  to  pinhole  geometry 


(a  is  the  pinhole  to  detector  distance,  b  is  the  pinhole  to  object  dis¬ 
tance)  indicates  that  the  pinhole  point-spread  function  (full-width  half¬ 
maximum)  at  the  imaging  volume  center  is  2.5  times  the  pinhole  diame¬ 
ter.  Figure  26.4  demonstrates  the  observed  and  expected  line  broadening 
due  to  the  non-magnifying  geometry. 

5-  Conclusions  and  Open  Issues 

Low  resolution  imaging  with  a  static  pinhole  geometry  and  energy- 
integrating  storage  phosphor  detectors  is  a  feasible,  low  cost  alternative 
for  high-throughput  evaluation  of  imaging  agents  and  for  determining 
optimal  timing  of  high  resolution,  small  animal  SPECT  studies.  High 
resolution  is  not  limited  by  the  intrinsic  resolution  or  sensitivity  of  the 
SP  detector,  but  is  due  to  limited  angular  sampling  with  small  numbers 
of  projections. 

The  current  prototype  possesses  detector  areas  that  allow  for  expan¬ 
sion  to  2  rings  of  90°  or  3  rings  of  60°  pinholes.  The  challenge  for  higher 
resolution  is  to  increase  angular  sampling  by  increasing  the  number  of 
projections  in  the  static  geometry.  Sixty  projections  has  been  reported 
as  sufficient  for  imaging  [Palmer  and  Wollmer,  1990],  but  multiplexed 
multi-pinholes  [Schramm  et  al.,  2001],  coded-apertures  [Sehellingerhout 
et  al.,  2002],  and  relaxed  collimation  [Prior  et  al.,  1993]  are  possible 
solutions  for  efficient  projection  packing. 

For  low  exposures  with  PS  image  plates,  it  has  been  shown  that 
noise  levels  obey  the  expected  Poisson  statistics  [Arbique  et  al.,  2003]. 
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However,  at  uniform  higher  exposures,  it  has  been  demonstrated  that 
structured  noise  source  in  the  image  plate  contributes  to  spatial- variant 
signal- to-noise  ratio  [Nishikawa  and  Yaffe,  1990].  A  procedure  for  cor¬ 
recting  projection  images  in  the  pre-processing  stage  has  not  been  ad¬ 
dressed. 

BaF(Br,I):Eu2'f  SP  imaging  plates  are  optimized  for  80  kVp  radio¬ 
graphy.  Unfortunately,  a  large  K-edge  at  40  keV  reduces  the  sensitiv¬ 
ity  of  SP  detectors  to  1251  emissions.  The  future  commercial  availabil¬ 
ity  of  storage  phosphors  more  suitably  matched  to  low  energy  emis¬ 
sions  (RbCsiTl4*  and  CsBr:Eu2+)  will  increase  the  efficiency  of  energy- 
integrated  imaging  with  125I. 

As  stated  above,  the  primary  challenge  remains  development  of  image 
reconstruction  algorithms  that  provide  scatter  correction.  Inevitably, 
this  will  probably  require  knowledge  of  the  mouse  topography  (size  and 
shape)  for  modelling  scatter  on  a  per-animal  basis. 
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Notes 

1.  The  terms  screen  and  image  plate  axe  synonymous  in  the  context  of  this  paper 

2.  The  terms  storage  phosphor  (SP),  photostimulablc  luminescence  (PSL),  and  photo- 
stimulable  storage  phosphor  (PSP)  are  synonymous  in  the  context  of  this  paper. 

3.  Exposure  time  will  depend  upon  collimation.  Our  experience  is  that  1  minute  per  2 
/iCi  is  typically  sufficient 
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